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Abstract 

We investigate the effect of a crystal edge dislocation on the metallic surface of a Topological 
Insulator. The edge dislocation gives rise to torsion which the electrons experience as a spin 
connection. As a result the electrons propagate along confined two dimensional regions and circular 
contours. Due to the edge dislocations the parity symmetry is violated resulting in a current 
measured by the in-plane component of the spin on the surface. The tunneling density of states for 
Burger vectors in the y direction is maximal along the x direction . The evidence of the enhanced 
tunneling density of states can be verified with the help of the scanning tunneling technique. 
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I Introduction 



Seldom a new state of matter is predicted. Even less often is it unambiguously observed 
in the laboratory in the spectacular way was done with HgTe quantum wells Indeed, 
this novel two-dimensional (2D) topological insulator (Tl) shows a quantized charge con- 
ductance, similar to that observed in Quantum Hall effect, without the need to apply a 
magnetic field. This result triggered further theoretical and experimental work that resulted 
in the discovery of metallic surface states in the insulating alloy Bi^Sbi-x, the first three- 
dimensional topological insulator. In this case, the conducting surface states are formed by 
topological effects that render the electrons traveling on such surfaces insensitive to scatter- 



ing by impurities. The theoretica 



of the Brillouin zone 



m 
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foundations of this phenomena are based on the topology 
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17, 



22 



3l|. The theoretical predictions have been con- 



firmed experimentally for the three dimensional {3D) Topological Insulators (TJ) Bi2Se5, 
Bi2Tes and Bii^^Sbx- The electronic band structure of these crystals is time reversal in- 
variant obeying Kramer's theorem = — 1 and have a single Dirac cone which lies in a 



gap 



23| . At the boundary of the 3D TI, one obtains a 2D surface with an odd 



number of chiral edge excitations coined helical liquid [37j and realized experimentally in 
the two-dimensional CdTe/HgTe/CdTe quantum wells [l|, [sgI. The dissipationless surface 
states are believed to be quantum-protected by the bulk insulator. A variety of transport 
experiments suggest that the conductivity of a 3D TI contains a significant metallic contri- 
bution 391]. Scanning tunneling microscopy {STM) and transmission electron spectroscopy 
{TEM ) of the chiral metallic boundary show that crystal defects modify the density of 
states. The momentum- resolved Landau spectroscopy of the Dirac surface state in Bi2Se^ 
reveals, in addition to the Landau spectroscopy, triangular-shaped structure caused by 



the presence of vacancies of Bi at the Se sites. These experiments suggest that the two 
dimensional surface of the TI is sensitive to defects and geometry and therefore, the quanti- 
zation rules are expected to be modified, and probably are revealed through the new Berry 
indices 7 ^ i |42|. (I n a recent Shubnikov-de Haas oscillation experiment performed on 
the TI Bi2Se3 42[ the authors show magnetic oscillations which correspond to the Landau 
level quantization En = Vp^y2{n + '-))heB with the Berry phase index which is different 



from 7 



2 V 



The appearance of the new indices 7 7^ 5 might be due to topology of the 



defects or/and interactions. Recently the effect of the 5*^^^(r) impurity potential on the 
metallic surface of the TI has been shown to modify the local density of states. The authors 
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in ref . 



19] have proposed that crystal dislocations generate protected zero modes which give 



rise to perfect metallic conduction. This idea has been used to explain the thermoelectric 
transport in the 3-D TI 25 1. 



From the theory of quantum crystal 



471] three primary types of topological defects are 



known : edge dislocations, screw dislocations and mixed dislocations. The strength of the 
dislocations is measured in units of of the Burger vectors which corresponds to the shortest 
lattice translation in the crystal. Experimentally, the dislocations are seen as dark lines in 
the lighter central regions of the transmission electron micrograph. 

The dislocations induce crysta 
shear modulus and Poisson's ratio 



stresses which are characterized by the Burger vector, 
3^, 47]. As a result, the lattice vibrations are complicated 



and the electron-phonon interaction plays a crucial role in inducing new dynamical effects. 
We will ignore the dynamical effects and concentrate on the static effects caused by the 
elastic strain field on the electrons. In a crystal we have the neutrality condition, namely 
the sum of the Burger vectors must be zero. We consider first a single dislocation (one 
assumes that a second dislocation with an opposite Burger vector is located far away from 
the first dislocation) and in the second stage we generalize the results to an even number of 
dislocations. 

Thirty two years ago the propagation of electrons in the presence of a dislocated crystal 
was investigated (26|. Those calculations, performed with the Schrodinger equation, revealed 
the possibility for interesting electronic transport. For a screw dislocation, the Schrodinger 
equation is equivalent to the Aharonov-Bohm problem in two dimensions [s^] which hints 
that the physics of persistent currents might play a crucial role . In graphene, different 
topological defects such as crystal disclinations have been modeled as a vortex field 



28 



29]. 



We start our discussion with the Peierls model [20| which is an approximation for the 
edge dislocations without torsion . This model gives results which are in agreement with 
those obtained by for the 3D TI. For this model we find that at the boundary of the 
dislocation a zero mode exist and confines the electrons to the dislocation line. 

Unfortunately, the Peierls model represents a over- simplification of the reality. The model 
is lacking torsion which is an important property for dislocations. A complete description of 
;he edge dislocation must include the torsion tensor . This tensor generates spin connections 
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44j-|46| which are controlled by the Burger vector. Using the complete description of 



the edge dislocation we find that the electronic excitations are confined to a two dimensional 
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region and to a set of circular contours. Such structure can be studied using the formulation 



introduced in ref. 50l-l52|. The circular contours have a radius Rg{n) n = ±1,±2,±3,... 
which are determined by the strength of the Burger vector. Comparing the results obtained 
from the model with torsion to the one obtained from the Peierls model, we find that the 
torsion destroys the zero mode state, but due to the Parity violation the system is only 
weakly affected by backscattering. 

The contents of this paper is as follows: In chapter II we present the chiral model for 
the boundary surface of the TI. Chapter III is devoted to the edge dislocation. We devote 
a short discussion to the Peierls model and present an extended derivation for the effect of 
dislocations on the TI. Chapter IV is devoted to the solution of the wave function in the 
presence of an edge dislocation. Chapter V presents the explicit wave functions for the two 
dimensional region and circular contours. Chapter VI presents the the computation of the 
tunneling density of states induced by the edge dislocation. In chapter VII we show that 
the violation of the parity symmetry by the edge dislocation generates a current which has 
an in-plane spin component. We suggest that this current which confirms the presence of 
the edge dislocation might be measured using Magnetic Force Microscopy techniques. 
Chapter VIII is devoted to conclusions. 

II-The chiral metal - the boundary surface of the three dimensional Topolog- 
ical Insulator 

The low energy Hamiltonian for the bulk 3D TI in the Bi2Se3 family was shown to 



behave on the boundary surface (the x,y- plane) as a two dimensional chiral metal 11 | 



which is similar the Rashba lOj model. 

(1) 

h^'\x, y) = hv F[i(y^ dy—ia"^ dx] is the chiral Dirac Hamiltonian in the first quantized language, 
fi? ~ 5 ■ 10^^ is the Fermi velocity, a is the Pauli matrix describing the electron spin and 
/i is the chemical potential measured relative to the Dirac F point. The Hamiltonian for 
the two dimensional surface L x L describes well the excitations smaller than the bulk 
gap of the 3D TI at 0.3 eV. Moving away from the F point, the Fermi velocity becomes 
momentum dependent; therefore, we will introduce a momentum cut off A to restrict the 



validity of the Dirac model. The chiral Dirac model in the Bloch representation takes the 
form: h = hvpiK x a) ■ z = hvpi—cr^ky + a^k^) where periodic boundary conditions imply 
= ^rn^ , = 0, ±1,±2, .. and ky = , niy = 0, ±1,±2, ... The eigen-spinors 

for this Hamiltonian are : \u{K) >= [\u^{K) >, \ui{K) = \K > ze^^^'^^'^'f)]^ where 
x{kx, ky) = tan~^{Y') is the spinor phase and e = hvp^^kl + ky is the eigenvalue for particles 
. For holes we have the eigenvalue e = —hvp^/kl + ky and eigenvectors \v{K) >= [\v^{K) > 
, \vi{K) >]'^ = \K > (g)[— 1, ie*-*^*^'^""'^''-']^. The chirality operator is defined in terms of the 
chiral phase x{kx,ky): 

K 

{a X — ^) ■ z = sm[x{kx, ky)](T^ - cos[x(fcx, ky)](T'^ (2) 
\K\ 

The chirality operator takes the eigenvalue — (counter-clockwise) for particles 
[sm{x{kx,ky))a^ - cos{xikx,ky))a'^]\K > ®[l,ie'^(''-'''yY = -\K > ®[l,ie'^^''-'''yY and 
+ (clockwise) for holes [sm{x{kx, ky))a^ - cos^xikx, ky))a^]\K > ze*'^^^^'^^'^]^ = \K > 

1, ■ie*^^''"^'''''^]^. The Hamiltonian h^-\K) = hvpi—a^ky + a'^k^) is time reversal in- 
variant h^'^'{—K) = TK^'^-{K)T^^ (T = {—iay)Ko is the time reversal operator and Kq 
is the conjugation operator). At = we have h^'^'{K = 0) = Th^'^'{K = 0)T~^ 
therefore the eigenstate is a Kramer degenerate. The eigenstate \u{K) > obeys 
T'^\u{K) >= —\u{K) >. The states \u{K) > and T\u{K) > are orthogonal to each other. 
Since T\u{K) >oc \u{—K) > the property = — 1 guarantees that backscattering is pro- 
hibited. 

III-The edge dislocation on the metallic Surface of the three dimensional 
Topological Insulator 

A-The Peierls model for the edge dislocation 



We start our presentation with the Peierls phenomenogical model [20| which might be 
useful for a large number of dislocations. An edge dislocation in two dimensions can be 
formulated in the following way: We cut the the two dimensional crystal into two halves 
{y < and y > 0). After cutting the crystal along the line y = 0, two surfaces [x,y < 0] 
and [x, y > 0] have been created. We glue back the two halves by translating by minus one 
half of the lattice spacing along x for ?/ < and one half along x for ?/ > 0. This procedure 
creates a dislocation along the line y = with a Burger vector along y , B^'^\ We introduce 



a distribution of dislocations p{x,y) then the quantity p{x,y)dxdy represents the Burger 
vector passing through the area dxdy. In order to understand the effect of the dislocations 
on the chiral fermions, we rewrite the two dimensional model on a lattice. On a lattice the 
model is characterized by the hopping matrix elements t = to (for opposite spins) and t' = tg 
( for parallel spins). In the absence of the Peierls deformation the lattice model takes the 
form: 

[x=na,y=ma] 

+ (-2t')^t(r)V^^(r) + (-2t')^t(r)Vy^(rO 

(3) 

where dy, represents the discrete lattice derivatives and the action of the discrete op- 
erators V^. Vy on the field operator is: Vx'^ir) = ^(r + a^) + \l/(r*— a^) — 2\E'(r) and 
Vy\E'(r) = \l/(r + tty) + \l/(r — ay) — 2\l/(r). When 2tQ = M we obtain the chiral model given 
in eq.(l). 

H^^ J2 'if^{f)[a''{-it)dy + ay{it)d^]^{f) (4) 

[x=na,y=ma] 



Following 20] we use the explicit form of the deformation field introduced by the dislocations 
u{y) = tan~^(X)). According to the elasticity theory, the Peierls edge dislocation [47J 
is determined by the half width W function and Burger vector B^"^^: p{y)dy = B'^'^\ p(y) = 
— u{y) oc M{y) = ^^^tan"^(^). Replacing the coordinates [x,y] by the deformed 

one, [X{r),Y{r)] we introduce the deformed matrix elements t' ^ : t' = + {-S^)u{y). 



As a result we find: -2t' + M = {-2t'n + M) - 2{p^)u{y) = - hvpnuiy) where hvp = 2to 

dt 

, K = j-fTj— As a result we obtain the Peierls Hamiltonian which is similar to the domain 



wall model introduced by |32|. 



e^{x,y) = h''-'-''\x,y)i,{x,y) 

h^^'^''^\x,y) = hvF[-i(y''dy + ia^d^ - a^KM{y)] 

(T^i'xi.x, y) = Xipxix, y); A = ±1; 

(5) 

The eigenfunction for the Peierls model obeys periodic boundary conditions in the x direction 
'ip{x + L,y) = ip^x^y). Therefore the eigen- spinor is given by : 'ip{x,y) = e*P^$(|/), p = 



^m, m — 0,±1,±2, .... The eigenvalue e is replaced by E where E = e — hvpp- We 
observe that if ip{x, y) is an eigenfunction with the eigenvalue E, the eigenfunction a'^il){x, y) 
corresponds to the eigenvalue —E. E = corresponds to the eigenvalue e = HvfP with the 
zero mode wave function ipo^Xjy^p) given by : 





kb(^) rv ^.,/j.„„-i/ «/' ^ / 1 

—i 



9[y] represents the step function which is one for y > and zero otherwise. As a result, 
we obtain for B^^^ > only a single zero mode with A = 1 for y < and A = —1 for 

y > 0. The eigenstate given in eq.(6), \Uo{p) > and \Uo{—p) > are related by the time 
reversal symmetry T\Uo{p) >oc \Uq{—p) >. From the equations = —1 and T\Uo{p) >oc 
\Uo{—p) > we conclude that the zero mode state is stable against backscattering along the x 
direction. ( The states \Uo{p) > and \Uo{—p) > are defined according to their wave functions 
<x,y\Uo(p) >='^o{x,y;p) , < x,y\Uo(-p) >= 'ijjo(x,y; -p).) 

(2) 

B-The edge dislocation with with a non zero torsion T^i^l 

For the remaining part, we will limit our discussion to the microscopic formulation of 
dislocations. When an edge dislocation is introduced into the crystal, the lattice coordinates 
f = {x,y) are modified, 'r^R = f+u = [X(f),Y{r)] where u{r) is the local lattice 
deformation with the core at the dislocation centered at 7^= (0, 0). The elastic field contour 
integral around the core f = (0,0) is determined by the Burger vector B'^"'\ a = 1,2: 
^ dx^dij,u^°'\f) — —B^°'\ For an edge dislocation in the x direction the Burger vector 
is in the y direction . The value of the burger vector B^'^^ is given by the shortest 
translation lattice vector in the y direction. (For the TI Bi2Se^ the length of the vector 
B^'^'> is 5 times the inter atomic distance ). For the two dimensional surface we use the 
notation .ji = x,y for the fixed Cartesian coordinate and X" ,a = 1,2 to describe 
the reference frame , the media with dislocations. In the media with dislocations we 
introduced a set of vectors Ca which are orthogonal to each other (e^, Ca) =< e''|ea >= S^. 
The unit vector Ca can be represented in terms of the Cartesian fixed frame space with the 
coordinate basis ,ii — x, y. The basis in the media frame can be expanded in terms of the 
fixed Cartesian frame d^; we have : = e^5^ (for the particular case where vectors are 



— * 

given by = da the transformation between the two basis is = (5^). Any vector X can be 
represented in terms of the unit vectors Ca (in the dislocation space) and in the Cartesian 
fixed coordinates space , X = X'^Ca = X^d^. The dual vector e" is a one form and can 
be expanded in terms of the one forms dx^. We have e"" = e^^dx, where represents the 
matrix transformation = {d^X°')dx^ . The scalar product of the components e^e" = g^^^^ 
e^^el = 5afi defines the metric tensors, g^^i, (in the Cartesian frame ) and 5^,6 in the dislocation 
frame. We will compute the matrix elements fields for our problem : 



el = d,X''{r); a = 1,2; ^i = x,y (7) 

Following 



33| we can express the Burger vector in terms of the the partial derivatives with 
respect the coordinates a = 1, 2 in the dislocation frame and fi = x,y for the fixed Cartesian 
frame : 

9.4 - S.el = B<-'^S\f) (8) 
Usiug Stokes theorem, we replaee the Hne integral j d.'eUr) by the surface integral 

j J dx'^dx'^ldxCy — dye'^]. For a system with zero curvature and non zero torsion T^fj we find 
that the surface torsion tensor integral / / dx^dx^Tj^} is equal to J J dx^dx"[dxe^ — dye"^], 
and therefore both integrals are equal to the Burger vector. 

d.el - d,el = Bm6'(f) 

(9) 

where dx^dx'^ represents the surface element. The tangent components can be expressed 
in terms of the Burger vector density B^'^^S'^{r) 3^ : 



2 ^^^^ y 2 B^"^^ ^ 



27r ' (x2 + y2) ' y ^ 2tx ' (x2 + 1/2) 



el = l- el = (10) 



Using the tangent components, we obtain the metric tensor g^ ^. 

^ elel + eje^ = ^^,,(f); e^^ ^ e^e^ + e^ej = 5,,, (11) 
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The inverse of the metric tensor g^^vi'f^ is the tensor g^'^{r) defined trough the equation 
gfi,T{r)g^''^{r) = S'^- Using the tangent vectors, we find to first order in the Burger vector 
the metric tensor (7^ and the Jacobian transformation \/G: 

_ B^'^ y B^^) y _ B^^) y 

g.,y- 2vrx2 + y2' 3y,y-^ 27rx2 + y2' 9y,.-yi. U - aet[g^,,\ - I ^2 ^ y2 

(12) 

The inverse tensor is given byigf^'^ ^ 1, g^'^ = g^'^ = — 4~ 2? 2 ? g^'^ = 1 + 2I 2 - Using 
the inverse tensor g'^''^ we obtain the inverse matrix which is given by: 

e^a = ea^ug"'' = {Kbet)g'''' = etg"'^ (13) 

In figure 1 we show the coordinate transformation f = {x,y) — )■ [X{r) = x,Y{r) = y + 
-^^tan^^(|)] with the core of the dislocation centered at f= (0,0). 

Next we consider the effect of the dislocation on the TI Hamiltonian. Using the compo- 
nents we compute the the transformed Pauli matrices. The Hamiltonian in the absence 
of the edge dislocation is given by h'^'^' = i'^°'da = X]a=i 2 ^7"^a where the Pauli matrices are 
given by 7^ = — ^ = cti and 7^ = cr^. (We will use the convention that when an index 
appears twice we perform a summation over this index.) In the presence of the edge disloca- 
tion, the term 7"9a must be expressed in terms of the Cartesian fixed coordinates fi = x,y, z. 
As a result, the spinor \l/(r) transforms accordingly to the SU{2) transformation . If \l/(-R) is 
the spinor for the deformed lattice, it can be related with the help of an SU (2) transformation 
to the spinor \E'(r) in the undeformed lattice: \1/(X, Y) = e~* '^V^'°'^vl>(x, y) . Where Sip{x, y) 
is the rotation angle between the two set of coordinates: 5(^{x,y) = tan^^(^) — tan^^{^). 
Using the relation between the coordinates (see figure 1) X = x, and Y = y + ^'tan'-^ {^) 
with the singularity aX x = y = gives us that the derivative of the phase which is a delta 
function, dxS(p{x,y) = —dy6ip{x,y) oc 6'^{x,y). Combining the transformation of the deriva- 
tive with the 5*0(2) rotation in the plane, we obtain the form of the chiral Dirac e quat ion 
in the Cartesian space (see Appendix A) given in terms of the spin connection cu^'^ |3l|: 

The Hamiltonian /i^'^' — t- h'^'^^'^ is transformed to the dislocation edge Hamiltonian with the 
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explicit form given by: 



fledge ^ ^^ig^ _ ^^2g^ ^ ta'ei^[d^ + l[a\ a'^f] - ta^e'ild^ + l[a\ a']cu'f] 



(15) 

To first order in the Burger vector we find : 00^^ = —u)"^ = and —(^y^ = oOy''^ = — -^^5^(r), 
see eqs. (69 — 71) in Appendix A. 

h""'' ^ ia\dy - ^a^B'^^^6\r)) - ta^d, (16) 

In the second quantized form the chiral Dirac Hamiltonian in the presence of an edge dislo- 
cations is given by : 

(17) 

j^edge jg Hamiltoniau in the first quantized language, /i is the chemical potential and 
\l'(r) = [\E'i-(r), \E'|(r)]"^ is the two component spinor field. 

IV- The solution for the metallic Surface in the presence of an edge dislocation 

(2) 

with the torsion T/^ j; 

A-The model in the momentum space 



We work in the momentum representation 12| where the edge Hamiltonian /z^'^^e takes 
the form: 

/i'^'^se ^ ;^^^[(^ xa)-z] + ^B^^^a^ j ^e*fe(*9^.)+5.(^9^.)l (18) 

The Hamiltonian is time reversal invariant but it is not invariant under the pla- 
nar parity symmetry. The eigenstate \u{K) > obeys the secular equation: [e — e{K) + 
i^a'^]\W{e) >= where \W{e) >= f -j^luiK - q) >. The the presence of the term 
generates unstable solutions . A stable solution can be obtained if the Pauli 
matrix annihilates the eigenstate. Using the chirality operator we observe that states 
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propagating in the x direction are eigenstates of o"^ (states polarized in the y direction). 
a^\K = K > (g)[l,ie*>^('=-'=^=°)]^ = sgn{K)\K = K > ie^>^('=-^^=o)]^. The eigen- 
states of the edge Hamiltonian which satisfy the chirahty operator must have the form : 
\u{kx,y) >= e^^'^^llu^ly) >, \ui{y) >]'^. In the next sections we will construct the explicit 
eigenstates for this Hamiltonian. 

B-Identification of the physical contours for the edge Hamiltonian h^^^'^ 



In order to identify the solutions, we will use the complex representation. The coordinates 
in the complex representation are given by, z = ^{x + iy) , z = |(x — iy), = — idy, 
9^ = dx + idy. In this re pre sentation the two dimensional delta function (5^(r) is given by 
5^(r) = ^dz{^) = 48|,|49(|. We will use the edge Hamiltonian /z^'^s^ and will compute 

the eigenfunctions Ue{z,z) = [Ue^{z,z),Uei{z,z)]'^ and Ve{z,z) = Ve^{z,z),Vei{z,z)]'^. The 



eigenvalue equation is given by: 



eU,^{z,z) = -[d. + (^)d,{^)]U,^{z,z) 

eU,,{z,z) = [d^+i^)d^i-)]U,,iz,z) 
yj In Z 



(19) 



The eigenfunctions Ue{z,z) and Ve{z,z) can be written with the help of a singular matrix 
M{z,z) y : 



Ue{z,z) = M{z,z)F^{z,z) 



g 2t: ^ z' 







13(2) .1. 
e 27r (z' 



F^{z,z) 
Fei{z,z) 

{F^{z,^) and F_^{z,^) are the transformed eigenfunctions for e > and e < respectively 
.) In terms of the transformed spinors the eigenvalue equation h'^'^^'^{z,l)u^{z,^) = eUe{z,l) 
and F^i{z,z) becomes: 





( F^iz,z)^ 




I{z,z) 




-d, 




{F^ii^,^) J 




o{nz,z)*_ 







where I{z, z) = 


e 27r \ zz 


= e vr y^-2+y-2) ^ ^^^^ ^-^y 


= e 




for zero modes e = and find : 

d,F,^{z,z) = 



c^F,^{z,z) = 



(20) 
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The solutions are given by the holomorphic representation F^=o^{z, z) = /^(-z) and the anti- 
holomorphic function F^=oi{z,z) = fi{z). The zero mode eigenfunctions are given by : 

«e=o,t(^) = e-'^^^^/tlz), u,=oM = e-'^^^^Uiz) (21) 

Due to the presence of the essential singularity at z = it is not possible to find 
analytic functions and fii'z) which vanish fast enough around z = such that 

J d"^ z{ue=o,\{z))*Ue=o,x{z) < oo. Therefore, we conclude that zero mode solution does 
not exists. The only way to remedy the problem is to allow for states with finite energy. 
In the next step we look for finite energy states. We perform a coordinate transformation 

z^W[z,z]; z^W[z,z] (22) 

We demand that the transformation is conformal and preserve the orientation. This restricts 



the transformations to holomorphic and anti- holomorphic functions |48|]. This means that we 
have the conditions dzW[z,z] = and dzW[z,z] = 0. As a result we obtain W[z,z] = W[z] 
and W[z,z] = W[z], which obey the eigenvalue equations: 

eF,^{W,W) = d^F,^{W,W) 

(23) 

This implies the conditions = {I{z,'z))* and ^^l^^ = I{z,^). Since I{z^'z) is neither 

holomorphic or anti-holomorphic and satisfy |/(2;,z)| = 1, the only solutions for W[z] and 
W[z] must obey I{z,z) = 1: 

J(z,z) = e^^^^^ = e*2™; n = 0, ±1, ±2.... (24) 

For I{z,^) 7^ 1 one obtains solutions which are unstable . The stable solutions will be given 
by a one parameter s curve (s is the length of the curve) r(s) = [x{s),y{s)] which obey the 
equation I{z,z) = 1. The curve r(s) allows us to define the tangent t{s) and the normal 



vectors N{s) SOj. This allows us to introduce a two- dimensional region in the vicinity of 



the contour of r(s) — )• R{s, u) = r(s) + uN{s). 
V- The wave function 
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A-The wave function for n — 

The condition I{z,z) — e ^^^+^' = 1 forn = is satisfied for y = and large value of 
y which obey 2^^(^2+p-) << 1 ■ The values of y which satisfy this conditions are restricted 
io I {z^z) = e ^^^+^' 1. This condition is satisfied for values of y in the range: 

2 (^^)<^<7<1 (25) 

We introduce the radius Rg = and find that the condition I{z,z) ^ 1 give rise to 
the equation for y. The solution is given by x'^ + (y ± ^Rg^ = (^RgY- Therefore, for 
\y\ > \d\ > {^)2Rg > 2Rg we have / pa 1 which corresponds to a free particle eigenvalue 
equations. 

€F^{x, y)=e ^ (-'+y'''> [-d^ + idy]F,^{x, y) 
~ + idy\F,i{x,y)\ 
eF,^{x, y)^e - [d, + idy]F,^{x, y) 

[d:, + idy\F^{x,y) 

(26) 

For \y\ > d the cigcnfunctions arc given by: Uc^-\{x.y) = e ^='+^y' Ff.^^{x,y), Ue,i{x,y) = 
e ^x-iy) p^^_^(^x,y) where F^^{x,y) and F^{x,y) are the eigenfunctions of equation (26). 
The envelope functions e ^tt '^x+iy-'^ g ^vr U-ij/-' which multiply the wave functions Fet(^) y) ? 
Fej^(a;, y) impose vanishing boundary conditions for the eigenfunctions U^,i{x, y) and U^^-]^{x, y) 
at y — >■ ±oo . therefore, we demand that the eigenfunctions U^^^{x, y), U^,i{x, y) should vanish 
at the boundaries y = ±-|. Since the multiplicative envelope functions for opposite spins is 
complex conjugate to each other we have to make the choice that one of the spin components 
vanishes at one side and the other component at the opposite side. Two possible choices 
can be made: 



C/,t(x,y = t)^e ^'^ WF,t(x,f) = C/a(x,y = -t)^e ^'^ ^^-(-^)'F^(x, ) = 



or 



b(2), 1 , s(2), 1 



U,,^{x,y = -f ) = e ^-+*(-^) V,t(x, -f ) = f/a(^>l/ = f ) = " ^^"^^ ^^(a:, f ) = 
Making the first choice, (both choices give the same eigenvalues and eigenfunction) we 
compute the eigenfunctions Fg]-{x, y) and Ff^{x, y) ioY \y\> d using the boundary conditions 
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Fe^{^. ?/ = ^) = 0; F,^(x, y = = (27) 

Due to the fact that the solutions are restricted to \y\ > d no conditions need to be im- 
posed at X = y = 0. In the present case we consider a situation with a single disloca- 
tion. This is justified for a dilute concentration of dislocations typically separated by a 
distance I ~ 10~^m. ( In principle we need at least two dislocations in order to satisfy 
the condition that the sum of the Burger vectors is zero.) The eigenvalues are given by 



e = ■^fivp^jP' + cp'. The value of p is determined by the periodic boundary condition in 
the X direction p(m) = = -^rn^ m = 0, 1, (iV — 2), (iV — 1) and a is the lattice 

constant a ~ X' "^^^ value of q will be obtained from the vanishing boundary conditions at 
y — The eigenfunctions F^^^{x^y) will be obtained using the linear combination of the 
spinors introduced in chapter II. In the Cartesian representation we can build four spinors 
Tp^q{x,y), Tp_g{x,y),T_p^q{x,y),r_p_g{x,y) which are eigenstates of the chirality operator 
and are given by: 



r,,,(x,y)=e^^-e 



rp,_g{x,y)^e'P^e- 



iqy 



iqy 



r_^,,(a;,|/) = e^^V^^ 



1 



1 




r_,,_,(a;,|/) = e-^'-e-^^ . _ (28) 



where tan[xip,q)] = J. 

The Hamiltonian h'^'^^'^ is not invariant under the symmetry operation : x ^ 
—x,a^ a^.a^ — > —a"^; [/i^^^"^, P^] 7^ therefore; we need to construct two independent 
eigenfunctions F^^^g^\x,y) for p > and f'^^^q^^ {x,y) p < 0. 

Fp'^^fhx, y) = A{q)Tp^g{x, y) + B{q)Tp^_q{x, y) 

y) = C{q)T_p,,{x, y) + L'(?)r_,,_,(x, y) 

(29) 
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Employing the boundary conditions given in equation (27) we obtain the amphtudes ^||y 
' ^Ify conditions for the momenta q. For the pair rp^q{x, y) , Tp_q{x, y) p > we 



obtain : 



g = g+ = yA; + } tan-i(^); A; = 1, 2, 3...; tan[x(p, q+)\ = (^) 

Li L/ p p 



(30) 



Similarly, for the second pair r_p ^(x, y),r_p _g(x, y), p > we obtain: 



F'^llp>Qq_)-^{x,y) = e-^^'^e-^^(^''«-)[e^(''-^+^^(^''«-^) + (-l)'=+ie-*(«-2'+^^(^''«-))]; \y\ > d 
Fei-P>o,q.U^^^y) = -ie-'^'^e-5x(p,9-)[e^(9-i/-^x(P,9-)) + (_i)fc+ie-i(9-y-ix(P,9-))]. > d 

q = q- ^ yk - Y tan-^(— ); A; = 1, 2, 3...; tan[x(p, ?_)] = (— ) 
L L p p 

e{-p, q-) = ±hvF\l {-Y^mY + q^ 



(31) 



For the state with zero momentum p = we find: 



^S=oS4(^' y) - ^2e^ cos[gy - ^]; |z/| > 
q = ^ + jk-k = Q,l,2,^... 
e{p = 0,q) = ±hvF\q\ 



(32) 

The eigenfunctions for the dislocation problem will be given ior \y\ > d by: ui"'~^'^\x, y) = 
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' ^ '^^ G^(a;,|/)L ^(-P>o,,-),t^^' 

t/r°'°'(.,s)»e-'e<A)f;tol(-.») 

(33) 

where G{x,y) = 1 — 1/2(3:2+^2) is the Jacobian introduced by the edge dislocation. The 
eigenstates are normahzed and obey: J dx J dy^yG{x, y){Ua"'~^'^\x, y))*ul7~^'^\x, y) 
5a,a'^ J dx j dy-yyG{x, y){Ua"'~^'^\x, y))*U^~^'^\x, y) t^^.a'- The normahzation factor 
2const.{Bi )) ^ ^ weak dependence on the Burger vector S(^) . This dependence is a 

consequence of the Jacobian which affects the normahzation constant (see appendix B). 

(The multiphcative factor e '^^ ^x±iv' gives rise to a weak non-orthogonahty between the 
states. This non-orthogonahty of the hnear independent eigenfunctions can be corrected 
with the help of the Grahm-Shmidt method.) 

For the present case, backscattering is allowed but it is much weaker in comparison to 
regular metals. This is seen as follows: Time reversal is not violated; due to the parity viola- 
tion, the eigenstates , y) ^uT ^'^\x, y) are not related by a time reversal symmetry 
{Tui"~^'^\x, y) 7^ u^e^~^'^\x, y)) . As a result, the backscattering potential Vp-p is controlled 
by a finite matrix element between states with different eigenvalues e(— p, g_) 7^ e{p, g+) (con- 
trary to regular metals where the impurity potential Vp-p connects states with the same en- 
ergy). In the present case |e(— p, g_) — e(p, g+)| = hvp] {^mY + q^_ — ^ (^m)^ + q^]] ^ 
the eigenvalues are not equal, therefore the finite matrix element controlled by the backscat- 
tering potential Vp-p gives rise only to a second order backscattering effect! 

B- The circuleir contours-the wave function for n 7^ 

_ 2 b(2) , iy(a) X 

The equation I{z,z) — e — e*^'^" gives the set of ring contours for n — 
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±1,±2, ±3, ... shown in figure 2. The radius Rg for the fundamental contour(n = 1) is 

lg = ^ and Rg{n) = 

{x{s)f + {y{s) ± Rg{n)f = {Rg{n)f (34) 



represented in terms of the Burger vector B^'^\ Rg = and R„{n) 



The centers of the contours are given by = [0,i?g(n)] for n 7^ 0. When n > the 

center of the contours has positive coordinates (upper contour) and for n < the center has 
negative coordinates (lower contour). Each contour is characterized by a circle with a radius 
Rg{n) = centered at [x = 0,|/ = Rg{n)\ (see figure 2). The contour is parametrized in 



terms of the arc length < s < 27r-j^ which is equivalent to < < 27r . Each contour i 
parametrized by r(s) = [x(s), Rg{n)+y{s)] where x{s) = Rg{n) cos[-^^^] = Rg{n) cos[ip\ and 
y{s) = Rg{n) sin[-^^^] = Rg{n) sm[(p\. We will extend this curve to a two dimensional strip 
with the coordinate u in the normal direction: For the curve curve r{s) = [x{s),y{s)] we will 
use the tangent t{s) and the normal vector N{s) 5^. Therefore, the two dimensional region 
in the vicinity of the one parameter curve r(s) is replaced by r(s) — > R{s, u) = r(s) +uN{s). 

s s 
X [s, u) = Rg{n) cos[ ] + ucos[ 



IS 



Rg{n) Rg{n) 

s s 
y{s,u) = Rg{n)sm[ ]+Msin[ 



Rg{n) Rg{n) 



(35) 



We will restrict the width \u\ such that e*^™e^*'' ~ 1 where t] obeys 77 < ^ < 1 , |m| < 
■^s^ _ Rg{n) ^ Rg{^){^^) < '^g^"'' ■ In these new coordinates, the Dirac equation is 

27m 

approximated for \u\ < Rgin){-;^) = by : 

+ R,{n) 
^ + Rain) 

(36) 



The solution for the contour n ^ 0, < s < 2TxRg{n)\ \u\ < 



Din) 



The periodicity in s allows us to represent the eigenfunctions in the form: F^^{s,u) = 
ET=-ooEge'^^^e^'^-F,^{j,q) and F,^(s,n) = ET=-ooE,e^'^'^^^^e^'^F,^{j,q). We 
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find: 



Kg[n) 
j + 1 

Kg[n) 



(37) 



The determinant of the two equations determines the relation between the eigenvalue e, the 
transverse momentum Q{e) and the eigenfunctions F^^{j,q),Ff^{j,q). The eigenvalues are 
degenerate and obey : e(j — l;k) — e{j — —{I -\-l);k) ,where I > 0. 



±g(e); Q{e) = Je 



I + i 
2 (__2_\2 



F,{1, q) ^ [F^{1, q), F,^{1, q)f oc [1, e-^'^W)]^. ^(g^ ^ ^^-^(^^ 



(38) 



The value of the transversal momentum Q(e) will be determined from the boundary con- 
ditions at ±-^^. We will introduce a polar angle 9 measured with respect the Cartesian 
axes: The angle < ip{n = 1) < 27r for the upper contour n = 1 centered at [a" = 0,^ = Rg] 
is described by the polar coordinate < ^ < tt measured from the center of the Cartesian 
coordinate [0, 0]. The lower contour centered at [x — 0,y — —Rg] characterized by the angle 
< (p{n — —1) < 27r is described by the polar angle 9 restricted to tt < ^ < 27r. We 
estabhsh the correspondence between (p{n — ±1) and 9: 

(p{n — 1) —29 -\ for the upper contour n — 1, < 9 <7r 

(f{n = —1) = 29 + — + TT for the lower contour n = —1, < ^ < tt 

(39) 

Following the discussion from the previous chapter we will introduce the following boundary 
conditions: 

F^-'\s^u^ ^) = 0; F^r'\'^V- -y) - 

F^=~'\s.u = = 0; F^--'\s,y= ^) = 

D(n = ±1) = D (40) 
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For the two contours n — ±1 we introduce eight spinors F^g ^^\ip{n — 

±1),«),f!!J^V(^ = ±1),«), rLtcf^H<^(^ = ±1),«), r57J^Q^(^(n = ±1),«). using this 
spinor we will compute the eigenfunctions: Fl^^^^'^\(p{n = l),u) , F^i~Q~^'^\(p{n = —l),u) 
characterized by the transverse momentum Q- and the pair i^/g" ^'^\(p{n = —l),u) , 
F^i~^'^\ip{n — l),u) characterized by the momentum Q+. 



r'i'^=^'\ip{n = ±l),u) = e^K¥'(n=±i))giQ« 



gj(^(n=±l))g-jK(Z,Q) 
1 

gi(<^(n=±l))giK(Z,Q) 
1 

_g-«(¥'("=±i))g«'^(^i<9) 



r^T^'\^(n = ±1), li) = e-^'(^("=±i))e-^«" ( ^ | (41) 

Using the vanishing boundary condition given in equation (40) , we find for the upper contour 
n — \ and positive angular momentum I the wave functions: 

(42) 

For the lower contour n — —1 and negative angular momentum we find: 

Fii-i'Q^Z^^i^ = = e-iMn=-i)el«{',Q-)[e^(Q-"-^('.Q-)) + ^_i^k+i^-i{Q-u-lK(i,Q.))-^ 

^i(-7,Q:u(<^(^ = -1).^) = e-'^^("=-^)e^'^(''«-)(-l)e-*^(-^)[e*('?-"+^'^(''^-)) 
+(_l)'=+ie~^('3-w+|«('.Q-))] 

(43) 

For both set of eigenfunctions we obtain the quantization conditions for the momentum Q- 
and eigenvalues e{l, Q-). 

Q- = ^k-^ tan-^(^^^), = 1, 2, 3.. .; tan[«:(/, Q.)] = (^y^) 



e{l,Q_) = ±hv^^l{^^r + Ql (44) 
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The presence of the singular transformation e ^x+iy' and e ^-x-ia-' demands that the 
eigenfunction should vanish at a; = y = 0, which corresponds to the point u = 0, ip{n = 
1) = — I for the upper contour and u = 0,(p{n = —1) = | for the lower contour. Using 
the mapping given eq.(39) we can transform from the variables (f{n — ±1) to the polar 
coordinate 9 (the wave function must vanishes at ^ = and 9 — tt). In order to obtain a 
finite wave function, we combine the spinors Fi^^'^\(p(n — l),u), F^i'^y^\(p{n — —l),u) 
with the singular transformation functions e ^tt '^x+iy> ^ g ^tt yx-iy>-_ 



U,^i,Q_)^{9,u) = e-^[(^«+")'=°^[^("=^)l+^(^«+(^«+")^^°[^("=^)ll">,|j:^f)((^(ri = 



and 

,(2) 



The amplitudes r(Z, t) are determined by demanding the vanishing of the wave func- 

tion Ue{i,Q_)^{9,u) and Ue{i,Q_)_\,{9,u) at x — y — 0. As a result, we obtain the explicit 
form of the wave functions U^(^i^q_)^{9,u) U^(i^q_)^{9,u) (see equation (76) in Appendix-C) 
which depend of the parameters rj{u) and C{9,u) : r]{u) = -^^f^^: < 1; C{(^:U) = 



27r(Rg(l)+u)((sin[26>])2+(,7(n)-cos[26l])2)- 

For the second pair FiQ^^'^\'p{'n = — , F^f~Q^\ip{n — we obtain the quan- 

tization conditions: Q+ = ^A: + ^ tan-^( ^+y k=l,2,3... tank(Z, g+)] = (^±^) with 

'+2 '"'"2 

the eigenvalues: e(/, Q+) = ±.hvF\J~{-^-^Y + Q+. Following the same procedure as we used 

for the pairs Fi^y^\(p{n = l),u) , F^i~^^'^\ip{n = —l),u) we obtain the wave functions 
Ue{i,Q^)-\{9,u) , U^(i^Q^)^{9,u) which are given in equation (77) Appendix-C. 

VI -Computation of the STM density of states 

The STM tunneling current / is a function of the bias voltage V which gives spatial 
and spectroscopic information about the electronic surface states. At zero temperature, 
the derivative of the current with respect the bias voltage V is given in term of the single 
particles eigenvalues: e(m, g_) = -^hvp^ (^m)^ -|- g^, e(m, g'+) = -^hvp^J {^mY -\- ,m = 
0, 1, 2, 3... for contour n = 0. For the upper and lower circular contours n — ±1, we have 
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:e(/,Q_) = ±hvF^{^Y + Ql ,e{l,Q+) = ±hvF^f{^y + Ql ,1 = 0,1,2,3... The 
STM density of states is computed for a voltage V between the STM tip and the sample. 
The tunneling current is a function of the bias voltage V and the chemical potential /i > 



30j: 



oc D{E = eV; s, m) = J] D^^'^E = eV; s, u) 



dV 



EE E Eif^'^"^°^"'''^^(^'^)i''^[^^+/^-^^^^\/(7^^)' + ^'] 



ri=± m qr=q+,q~ 



+ EE E T.\u!r=^''^'''Ho,^^mev+^,-vhvJ^ 

n=±l I Qr=Q+,Q- a V 

(45) 

(77 = + corresponds to electrons with energy < e < /i and rj = — corresponds to electrons 
below the Dirac point e < 0. For the rest of this paper we will take the chemical potentials 
to be /i = 120m (this is typical value for the TI ). We will neglect the states with t] = — 
which correspond to particles below the Dirac cone ( see Appendix -B). The density of states 
at the tunneling energy eV is weighted by the probability density of the STM tip at position 
[x,y] for n=0. The contours for n = ±1 will be parametrized in terms of the polar angle 
6 and transverse coordinate u. The proportionality factor J for the tunneling probability 
(not shown in the equation ) = JD{V; x, y) is a function of the distance between the tip 
and the sample. The notation D'^"^(V"; x, y) represents the tunneling density for the different 
contours. 

A- The tunneling density of states for the Peierls model D^'^'^'^'^'-'^iV; x,y) 

We consider first the Peierls model introduced in section /// — A. This model has a 
zero mode for which we have computed the wave function in equation (6). We find that the 
tunneling density of states density is confined to the quantum strip y = with a varying 
width W, \y\ < W determined by the dislocations distribution 

W W hvp L 

(46) 

In figure 3 we show the tunneling density of states for the parameters ^^^^^ = 0.1 in units 
of the half width W. 
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B-The tunneling density of states (V; x, y) for n — Q 



Summing up the single particle states weighted with occupation probability 
we obtain a space dependent density of states for the two dimensional 
boundary surface < x < ^ and the coordinate y is restricted to the regions | < < f 
and ^ < y < |. We will perform the computation at the thermodynamic limit, namely 
we replace the discrete momentum j^k by Y = and by X = ^ where N = ^. We 
find for the dimensionless momentum q = qa the equations : q±{Y) — nY ± tan~^[^^] 
where 27rX — pa — p. As a result we obtain the following density of states ^ 

M+._i ^ 1 ql - jjp 
^dY^ TT ql+p'' 

^dY^ TT f_+p'' 

(47) 

Using this results, we compute the tunneling density of states in terms of the energy /i + eV 
measured with respect the chemical potential and the transverse energy e±_ = hvFq±- 

+ y - -_] - + eV^ - E^ax]) ■ ((cos[ ^^^^^ ^ - + (cos[ ^^^^^ ^ 1/ - 
-) ( ) — p - 1 ,2+^2+,. J |- / VP- ; 



^/i^^F L 4:^G{x,y) Jo v^(/i + eV) 



e 



2 



+^(^[, + y _ ^] _ + - ^_]) . ((cos[i^. - J)^ + (cos[^^. - J)^)] = 



^hvp L 4^G{x,y) L ^ 2L 

for \y\ > d 



H[n + eV — is the step function which is one for fj, + eV — > and zero otherwise, 
a ^ is the short distance cut-off and E^ax = hvpA < O.SeV is the maximal energy which 
restricts the validity of the Dirac model. We observe in the second line that the asymmetry 
in the density of states 1 ± \ lI^^^v) '\J^ ~ ^^^^ ) cancels. 
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(48) 



Equation (48) shows that the tunnehng density of states is hnear in the energy fx + eV (in 
the present case we have looked only for energies above the Dirac cone ). For the chemical 
potential /j, = 120mV, the zero energy corresponds to the Voltage V = —120mV. The 
tunneling density of states has a constant part at energies ^jf- ~ 0.2mV for —120mV < 
V < —lld.SmV . For V > —lld.SmV the density of states is proportional to /i + eV. 

In figure 4 we have plotted the tunneling density of states as a function of the coordinates 
X and y. The shape of the plot is governed by the the multiplicative factor e " ^x±iy> which 
governs the solutions in eq.(35). We observe that the density of state is maximal in the region 
\y\ < 10^(2). 

Figure 5 shows the dependence on the voltage V and coordinate y. Wc observe the linear 
increase in the tunneling density of states which is maximal in the region \y\ < 10B^^\ 

C-The tunneling density of states D'^"'^^\V,x,y;fi, ..f2M) for 2M dislocations. 

For many dislocations which satisfy Yl'^i -B*-^'"'-' = ( sum of the Burger vectors is zero 
) with the core centered at [xu,,yw] ,w = 1,2..2M the coordinate r = {x,y) — )■ [X{f),Y{r)] 
is replaced by [X{r) — x,Y{f) — y + ^^^^ tan~^(^5f^)]- Following the method used 
previously, we find the edge Hamiltonian with many dislocations takes the form: 

. 2M 

h^'^^^w = 1, 2...2M) ^ ia^[dy a^B^^''"^5\f- fj] - ia^d^ (49) 

2 

w=l 

As a result, the wave functions are given by: 

^(n=0,^=l,2...2M)^^^y^ OC J] e'^^ (^-™-^)F^("=°) (x, t/) 

w=l,2...2M 

jjin=0,n,=l,2...2M)^^^y^^ "Q )^^(n=0) (53) 

w=l,2...2M 

Using these wave functions, we find that the tunneling density of states is given by: 

L'("=°)(V,X,y;fi,..f2M) OC [[ e - ^(.x-x^y'+(y-y^y'+a^> (51) 

M)=1,2...2M 

In figure 6 we show the tunneling density of states for an even number of dislocations in the 
y directions which have the core on the y — axes {fu, — [x^, y^, — 0], w — 1,2, 3, ...2M). 
We observe that the tunneling density of states is confined to the axes y — and resembles 
the structure obtained from the Peierls model given in figure 3. 
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D-The tunneling density of states D^'^-^^\V;9,u) for the n — ±l contours. 

Following the same procedure as used for the n — case and the eigenfunction given in 
Appendix-C, we find for the tunneling density of states: 

For the even /c's, we solve for the momentum (5+ and (5_ and find: 



f R(2)^2 °o / / _|_ 1 

27iRg{l)D{l)^yG{9,u) V ^9^^) 

2(-l)' sm[QrU + Qr)] sin[QrM - ^/«(/, Qr)] ■ 

(cos[Z(e + — ) - C(^, n){-r]{u) + cos[2e])] - cos[(Z + l)(e + — ) + C(e, u){-t]{u) + cos[2e])]]; 

(53) 

Similarly for the odd /c's we find: 



27ri?,(i)i)(i)yG(M)Q.4;,Q_t^ V'^.(i)' 

[(g-2C(Msin[2e] ^ g2C(Msin[2e])^^Qg[g^^ _ g^)])2 ^ (cOs[g,« + \k{1, Q,)\f) + 

2( — 1)' C0s[(5r1i + ^'«(^) Qr)] C0s[(5j.li — ■^k(/, Qr)] " 

(cos[/(^ + — ) - C(^, n){-r]{u) + cos[2^])] - cos[(/ + 1)(^ + — ) + C(^, u){-r]{u) + cos[2^])]] 

(54) 

For the present case the energy scale of the excitations is governed by the radius -Rg(l) and 
width D. The spectrum is discrete and we can't replace it by a continuum density of states 
as we did for the case n = 0. 

In figure 7 we show the tunnefing density of states at a fixed polar angle ^ = | as a 
function of the voltage V . We observe that the density of states is dominated by high 
energy eigenvalues. This solutions are localized in energy. The range of the spectrum is 
above [i + eV > 200mV which is well separated from the low energy spectrum controlled by 
the n — contour (which ranges from —120'mV to TOmV). 
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Figure 8 shows the tunnehng density of states as a function of the polar angle 6 for a 
fixed energy . The periodicity in 9 is controlled by the discrete energy eigenvalues. 

In figure 9 we show the tunneling density of states at a fixed voltage y as a function of 
the polar angle < 9 < tt and width \u\ < 0.1. 

VII-The charge current-the in plane spin on the surface 

A-The current in the absence of the edge dislocation 

Prom the Hamiltonian given in equation 1 we compute the equation of motion for the 
velocity operator: ^ = j^^[x,h] = Vpa^ , ^ = jj^[y,h] = —Vpa^. We multiply the velocity 
operator by the charge (— e) and identify the charge current operators : Jx = {—e)vF<y^, 
Jy — {—e){—VF)cr^- This also represent the '"real"' spin on the surface. Therefore, the 
charge current is a measure of the in-plane spin on the surface. 

Integrating over the y coordinate we obtain the current /J"^" in the x direction. Using 
the eigenstates rp^g{x,y) and T_p^g{x,y) of the h^-^- Hamiltonian 



we find {Tp^g{x,y)){a^){Tp^g{x,y)) = -{T_p^g{x,y)){a'^){T_p^g{x,y)) therefore, we conclude 
that the current /J-^- = is zero. 

B-The current in the presence of the edge dislocation 

We will compute the current in the presence of the edge dislocation. The current operator 



quantization form for the current density. The operator is defined with respect the to 





jedge^^^ y) will bc givcu in terms of the transformed currents. We find that the current 
density operator J^'^^^{x,y) is given by: 
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shifted ground state |/i >= |0 > with the energy E — e — /i measured with respect the 
chemical potential and spinor field 'i>n=oix,y)- 

Jf^x.y) =< ^x\<^>i^,{x,y)Jl^^'{x,y)<^!n=o{xM^i > (56) 

Using the spinor eigenfunction given in equation (35) and the second quantized form given 
in Appendix -B we find : 



E>0 



The current is a sum of two terms computed with the eigen spinor obtained 
in equation (33): [C/f =°'''^(x, y), C/f =°'''^(x, y)]^a2[C/{"=°'^)(x, y), C/f and 
\p^^~^'^\x,y),U^^~^'^\x,y)Ya^\U^~^'^\x,y),U^^~^'^\x,y which have opposite signs. 
Due to the parity violation caused by the dislocation, the density of states is asymmetric 
1 ± ^ L(p+v) ~ ^J^tv^"^ ) resulting in a finite current. We integrate over the transversal 
direction y and obtain the edge current /^=0'e<ige 

^— j ^ dy<i^\r''<^^{x,y)\i^>^ 



2 



-S(2) 



—e)vF L ^2/ 1 \ dye " ^x^+y'^+a^' r r r — - — {hvpjV) ■ e 



An 'hvF^ 'L\J_L L J_L L ^/G{^) J "J " V ^ ^ ^(q|)^ + (e±) 



\2 



H[ijl — ^/(eii)^ + {^±y ] is the step function which is one for ^/(eii)^ + (e±)^ < A*- The 
single particle energies are e± — hvFq± and e|| — hvpp. For L 10~^m, chemical potential 
/X = 120mV and ~ 100 we find that the current /_j^'=o,e'^s'e -g -^-^ ^-j-^^ range of mA. 

To conclude, we have shown that the presence of an edge dislocation gives rise to a 
non-zero current which is a manifestation of the in-plane component of the spin on the two 
dimensional surface . Therefore a nonzero value 7^=0'e(ifle _^ q indication of the 

presence of the edge dislocation. This effect might be measured using a coated tip with 
magnetic material used by the technique of Magnetic Force Microscopy. 
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(58) 



Vlll-Conclusions 



Using the full description of the edge dislocation in terms of the torsion tensor, we have 
shown that the singularity at the core center eliminates the zero mode . As a results only 
weak backscattering effect is allowed . Using this formulation we have shown that for the 
case n = the tunneling density of states is confined along y = 0. At high energies the 
tunneling density of states is confined to circular contours governed by the Burger vector . 
For a large number of dislocations we obtain a result similar to the one obtained from the 
Peierls domain wall model. The in plane spin orientation is a manifestation of the parity 
violation induced by the edge dislocation. We propose that scanning tunneling and Magnetic 
Force Microscopy are advanced experimental techniques which can verify our predictions. 
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Appendix -A 



We consider that a two dimensional manifold with a mapping from the curved space X", 
a = 1,2, to the local flat space x^, fi = x,y exists. We introduce the tangent vector [45| 
^/l(^) — ^^Jm^' = which satisfies the orthonormality relation e^{x)e^^{x) = 5a,b (here 
we use the convention that we sum over indices which appear twice). The metric tensor for 
the curved space is given in terms of the fiat metric 5a,b and the scalar product of the tangent 
vectors: e^(a;)ej!;(a;) = gn^i,{x). The linear connection is determined by the Christoffel tensor 



"^aA = -r;,.5A (59) 
The Christoffel tensor is constructed from the metric tensor gfj„u{x). 

^i-^ = ~\Y1 9^'^(^)l9'^9uA^) + df,guA^) - drQp^A^)] (60) 

T=x,y 

Next, we introduce the vector field V = V^da = V^d^ where a = 1, 2 are the components 
in the curved space and fi = x,y represents the coordinate in the fixed cartesian frame. The 
covariant derivative of the vector field V"' is determined by the spin connection o;^^ which 
needs to be computed: 



D^Vix) 



d,V''ix)+Lo^^,V'' 



(61) 



For a two component spinor, we can identify the spin connection in the following way: 
The spinor in the the curved space (generated by the dislocation) is represented by ^E'(X) 
and in the Cartesian space it is given by is given by \l'(x) 53|. The two component spinor 
represents a chiral fermion which transform under spatial rotation as spin half fermion: 



^(X) = e^'^^-2'"«*(f) 

(62) 
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We have used the anti symmetric property of the rotation matrix uja,b = —^b,a, and the 
representation of the generator S'^"'' in terms of the Pauh matrices. 
Therefore for a two component spinor we obtain the connection: 

D.^ix) = id, + ^a;^-''S,,,)vI/(f) ^ (d, + ^c^^'^ia,, a,])^ix) (63) 

Next we will compute the spin connection u'^'^ using the Christoffel tensor. In the phys- 
ical coordinate basis the covariant derivative D,V^{x) is determined by the Christoffel 
tensor: 



D,V''{x) = d,V''{x) + V'^^^V' (64) 

The relation between the spin connection and the linear connection can be obtained from 
the fact that the two covariant derivative of the vector V are equivalent. 

D^V' = elD^V^ (65) 
Since we have the relation = e'^V" it follows from the last equation 

D,[el] = D^d^e'' = {D,d,)e^ + d,{D,e'') = (66) 

Using the definition of the Christoffel index and the differential geometry relation Vg^di, = 
— r^^^c^A [45|], we obtain the relation between the spin connection and the linear connection: 

D,K] = d.elix) - ri,,e^(x) + co^^^eUx) ^ (67) 

Solving this equation, we obtain the spin connection given in terms of the Burger vector. 
We multiply from left equation by the tangent vector e^, replace F^^, by equation (54) use 



the metric tensor relations e°(x)e|^(x) = 6a^b, e?^(^)e"(x) = g^^^{x). As a result, we find j45| : 



■\eP'''e'''\dpe„,,-d^ep,,)el (68) 



We notice the asymmetry between e'^''* and Cay. 
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For our case we have a two component the spin connection u;^^ and oUy"^ 



(69) 



These equations are further simphfied with the help of equations (11 — 13) with = , 
— \ and the Burger tensor d^e^y — dycl — B^'^^S'^{r) . 

\B^^)mr)\-^-^ - (1 - f^^)^)(^^)(l - ^ 



is(^)5(^)(r-)[- 



5(2) 2y - 



X . 



2 - 27r 1/2 + a;2' 

(70) 

and 

= ^e'^'^(9,e2 - d^e\) - - d^e\) - \e^^'e^\d,{b^^,e\) - a.(5,,,ej)]e^ = 

\9''^^A{^^l - d.el] - Ig'^'^elidyel - d^y] - \g'''''elg'^^^el[d,el - d^e^Yy = 

R{2) d(2) ry{2) 

-^5(2)(f)^^'M - ^^^'Hrlb^V^ - ^^'Vleje2e; ^ -^5^'\r) 

(71) 

To first order first the Burger vector 5^2) the spin connections are given by : io^ = —io'^ ~ 
and ujf = -ujf ^ -\B^^^5\r). 

Appendix-B 

The spinor field operator '^[z,^\t) is decomposed into a sum with different localization 
contours n, ^(2;,z;t) = X^„=o±i '^n{,z,z]t).Th.e energy levels for n 7^ are controlled by 
the the inverse of the Burger vector . 

We will consider the case n = 0. For = 0, we have the expansion in terms of the 
eigenspinor u^{s,u) for electrons with positive chirality v^{s,u) for holes with negative 
chirality. 
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ds J du^J G{s, u)ul{s, u)ue'{s, u) — S{e, e') 
ds J du\/ G{s, u)vl{s, u)v^/{s, u) = S{e,e) 

OO 6<0 V<^(«'«) 

(72) 

Using the eigenfunctions u^{s,u), v^{s,u), we construct the field operator 'i/n=o{s,u;t) as 
a superposing of particles and holes. 

e>0 e<0 e>0 

(73) 

where is the annihilation operator for particles with energy e and feL^ is the creation 
operator for a hole with energy — e. and b-e annihilates the ground state |0 > ,Ce|0 >= 
b-e\0 >= 0. The operators obey anti-commutation relations: [c£,c^,]+ = = 
5{e,e'). 

The material properties of the topological insulators are such that the Fermi energy /i is 
positive . As a result we have electrons and holes with positive chirality, and deep holes with 
negative spin chiralities. The energy is measured with respect to the chemical potential /i, 
E = € — iJ, .As a result one obtains a shifted ground state l/x >= |0 > . 

We define new operators using the holes operators b^.^ and electron operator c^.^: 
CE+t,-+ = as for E > 0, ce+^-+ = /^l^; for -/x < E < and 6L(^_^^)._ = for 
— 2/i < E — /I. These operators annihilate the ground state : q;e|0 >= 0, /3_e|0 >= and 
7_B|d >= 0. 
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E=Ea E=0 -fJ, 

E=0 E=-n -H-Ea 

E=Ea E=0 E=-/i 

= Y aEUE+^,{s,u)e-'^* + ^ (3-EUE+t,{s,u)e-'^^ + ^ -flEVE+^,{s,u)e-'^^ 

E=0 E=-n E=-^i-Ea 

E=Ea E=h E=ij, 

= aEUE+n{s, ^)e~' fi* + PI;U-e+h{s, u)e'~f^* + ^ 7^v_£+^(s, ii)e' ft* 

E=0 E=0 E=fi+EA 

(74) 

where -Ea is energy the below the Dirac point . As a rcsuh, the edge Hamiltonian is given 
by: H^'9e ^ J=^^A ^t^^^E + PII^eE + E^J+s, ^^eE] 

For most of the cases, the chemical potential is large, and we can approximate the spinor 
operator by a sum of states for particles and holes, ignoring the deep hole band: 

*„=o(s, u; t) ^[Q;£iiE+^(s, u)e~'^^ + PI;U_e+^(s, u)e'^ (75) 

B>0 
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Appendix-C 

The wave functions are given by: 

+ (_l)'e-C(Msin[2e]g-<{e,«){-r,H+cos[2e])g-a(2e+^) ^^^^Q_y^ _ g_)]]. 

+ (_l)'e-^(^'")-°P^]e-^^(^''')(-^('^)+^°^P^l)e-^'('^+^) cos[g_M - Q_)]]; 

_(_l)'e-^(^.-)-"P'^]e^^('^'«)(-^W+-«[2e])g-^((^+l)(2e+f )) ^ i^^^^ 

_(_l)/e-C(^.«)--[2^]e'C(^.")(-'?(")+-«[2e])g-i((Hi)(2e+^)) cos[Q_ii + g_)]]; 

(76) 

Similarly for the second pair we obtain the wave function: 

^gC{e,«)sin[2e]g-^C(e,«)(,H-cos[2e])g-^^(2e+^) ^-^^q^^ _ 1 ^(^^ g^)]]. 
^gC{e,«)sin[2e]g-^C(e,«)(r,M-cos[2e])g-^^(2e+^) cos[Q+« - ^k{1, Q+)]]; 

U^l^Q^f^ie, u) - 2ie^'^(''«+) [_(_;L)^e-^(^'")«'"^P^le^^(^'")(-^(")+^°^P^l)e^('+^)(2^+^) sin[g+ii - -k{1, g+)] 

^gC(e,n)sin[2e]g^C(^,»)(r,W-cos[2(?])g-^a+l)(2(?+f ) ^ 1^^^^ g^^]]. 

2 

^SS!u(^,«) = 2e^'^(''^+)[-(-l)^e-^(^'")^'°P^le<(^'")(-''(")+'=°^P^l)e^^ - -k{1,Q+)] 

' 2 

^gC(e,n)sin[2e]giC(e,«)W«)-cos[2e])g-i(;+l)(2e+^) cos[Q+ti + g+)]]; 

2 

(77) 
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where G ^ (9, u) is the Jacobian transformation induced by the metric tensor. 
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35 



y(x 



)/B 




FIG. 2: The contours {x{s)f + {y{s)-^)'^ = {^f for n = ±1, ±2, ±3(in decreasing size ),i?g(n) = 
n = corresponds to the equation y{s) = and |y| > d (see the text). The the distance is 
measured in units of the Burger vector B^'^\ 
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-100 
100 



FIG. 3: The tunenline density of states for /)-Pe«eris/^ Peierls dislocation confined to 

° W W 

the hne |x| < W for the STM voltage V, fi + eV <0 
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FIG. 4: The tunneling density of states forn = , ^ cx Z)("=°)(;g^, fi = UOmV). The right 
corner represents the intersection of the x coordinate which runs from 30 (right corner) to —30 and 
the y coordinate which runs from —30 (right corner) to 30 in units of the Burger vector. 
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D(x/B=-2;y/B,V;//=120) 



20 



FIG. 5: The tunneling density of states for n = as a function of y and V ^ D^"" ^\^2) — 
—2, -^(2)] ^ = 120mV). The voltage range is —120 < y < 50 and the y coordinate is in the range 
-30 < ^ < 30. 
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FIG. 6: Many Dislocations - with the core of the dislocations at [xw,y = 0] , w = 1,2...2M; 
The maximum of the tunneling density of states is confined along y = 0. The coordinates of the 
tunneling density of states are restricted to : —40 < < 40 and —20 < < 20. 
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D(6'=Pi/2,u/B=0.01,)U=120 ;V) 
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FIG. 7: The discrete tunneling density of states for n = 1, as a function of the voltage V 
Di-=^){V;e = ^,^,fi = 120mV) 
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D(6';u/B=0.01,V=280) 




FIG. 8: The tunneling density of states as a function of 6 D^'^=^\e] ^ = 0.01, V = 280mV, fi = 
120mV) 



42 




43 



M.Konig et al., Science 318,766 (2007) 

B. A.Volkov and O.A. Pankratov, JETP LETT, vol.42, 179(1985) 
F.D.M. Haldane Phys.Rev.Lett.61, 2015(1988). 

Karl Jansen, hep-lat/921203. 
David B.Kaplan liep-lat/920601. 

Michael Creutz and Ivan Horwath Phys.Rev.SO, 2297(1994) 

C. L. Kane and E.J. Mele Phys.Rev. Lett. 95 226801 (2005) 
C.L. Kane and E.J. Mele ,'Pliys.Rev.Lett. 95,146802(2005). 



Liang Fu and C.L. Kane cond-mat/0611341 



Bychkov, Y.A. and Rashba,E.L J.Pliys.C17,6039-6045 (1984). 
H.Zhang et al. nature physics 5,438,(2009). 
D.Schmeltzer,Phys.Rev.B 73,165301(2006). 



J.E.Moore and L.Balents, cond-mat/0607314 



Andrew M.Essin and J. E. More, cond-mat/0705.0172. 

Xiao-Liang Qi, Taylor Hughes and Shou-Cheng Zhang , Phys.Rev.B78, 195424(2008) 
Xiao-Liang Qi and Shou-Cheng Zhang cond-mat/1008.2026 

Andreas P. Schnyder, Shinsei Ryu, Akira Furusaki, Andreas W.W. Ludwig cond- 
mat/0803.2786. 

M.Z.Hasan and C.L. Kane, cond-mat/1002.3895 

Ying Ran et.al , Nature Physics vol5,298 (2009) 

R.E. Peierls ,Proc.Roy.Soc. 52,34(1940) 

Chao-Xing Liu et al. Physical Review B 82 045122 (2010 "' 

D.Schmeltzer , In Dark Energy :Theories, Developments, and Implications ISBN 978-1-61668- 
271-2 "'Topological Insulators- Transport In Curved Space"' (chapter 10), Editors:K.Lefebvre 
and R.Garcia, pp. 1-24, Nova Science Publishers, Inc (2011). 

Y.Xia ,D.Qian, D.Hsieh, L.Wray, A. Pal, A.Bansil, D. Grauer, Y.S. Hor,R.J.Cava and M.Z. 

Hasan ,cond-mat/ 0908.3513. 

D.S. Novikov Phys.Rev.B. 76 ,245435 (2007) 

O.A.Tretiakov et al Cond. Mat/1007.2966 . 



44 



[26] K.Kawamura Z.Physik B 29,100-106 (1978) 
[27] Y.Aharonov and D.Bohm Phys.Rev. 115 (1959). 

[28] F.de Juan,A.Cortijo and M.A.H.Vozmediano, Phys.rev.B 76 165409(2007). 

[29] A.Cortijo and M.H.A. Vozmediano, Nuclear Physics B 763[FS](2007) 293-308. 

[30] C.Kittel , Introduction to Solid State Physics ,eight edition 2005 John Willey and Sons,Inc. 

see pages 524 — 525 and 528 — 531. 
[31] M.Nakahara '"Geometry, Topology And Physics "' Taylor Francis Press 2003. 
[32] R.Jackiw and J.R. SchriefFer,Nuclear Physics B190,253-265,(1981) 
[33] H.S. Seung and D. R. Nelson Phys.Rcv.B 38,1005 (1988) 

[34] H.Kleinert " 'Multivalued Fields in Condensed Mattewr ,Electromagnetism, and Gravitation" ' 

World Scientific (2008) pages 348-350. 
[35] Z.F.Ezawa," 'Quantum Hall effects'" World Scientific (2008). 
[36] B.Andrea Bernevig Taylor L.Hughes, and Shou-Cheng Zhang cond-mat/061139. 
[37] Wu.C.B., A. Bernewig and Shou-Cheng Zhang Phys.Rev.Lett. 96 106401 (2006). 
[38] T.Hanaguri et al. Phys.Rev.Lett. 82, 081305 (2010). 
[39] A.A.Taskin and Yoichi Ando , Phys.Rev.B 80,085303 (2009) 

[40] Peng Cheng, Canli Song, Tong Zhang, Yanyi Zhang , Yilin Wang, Jin-Feng Jia, Jing 

Wang,Yayu Wang , Physical Review Letters 105, 076801 (2010) 
[41] H.Peng et.al. Nature Materials ,9 225 -229 (2010) 
[42] J.G.Analytis et al. cond/mat 1003.1713. 
[43] R. Biswas and A.V. Balatsky, Phys.Rev.B 81, 23405 (2010). 
[44] Ayelet Pnueh 'J.Phys. A: Math.Gen.27 1345-1352 (1994) 

[45] M.B.Grccn,J.H. Schwartz and E.Witten, "'Some Differential Geometry" 'Superstrings Theory 
vol.2 Cambridge Monographs on Mathematical Physics pages 271 — 277 (1987). 

[46] N.D.Birrell and P.C. Davies '"Quantum Fields in Curved Space'" pages 81 — 85 Cambridge 
University Press (1982). 

[47] Arnold M. Kosevici '"The Quantum Crystal"' Wiley -VCH (2005). 

[48] P. Di. Francesco. P. Mathieu and D. Senechal "'Conformal Field Theory"' page 119, Springer 

Text in Contemporary Physics (1997). 
[49] The Holomorphic representation of the delta function was brought to my attention by colleague 

V.P.Nair. 



45 



[50] P.Exner and P.Seba J.Math. Phys.30 2574-2580 (1989) 

[51] J.Goldstone and R.L. JafFe, Physical Review B45,14100(1992) 

[52] R.C.T. da Costa,Physical Review A23,1982(1981) 

[53] M.Maggiore, "A Modern Introduction to Quantum Field Theory "' Oxford University Press 
(2005), pages 26-31. 



46 



